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Abstract 

. I put forward a model with vectorlike isosinglet quarks in which soft 

\q [ breaking of CP leads to strong CP violation only arising at two-loop level. 

CN 

Non-perturbative effects in Quantum Chrome-dynamics (QCD) may lead to 
P and CP violation, characterized by a parameter 9, in strong interactions. The 
experimental upper bound on the electric dipole moment of the neutron requires 
9<2 x 10~ 10 10. (Electric dipole moments of elementary particles are both P- 
and T- violating; in normal field theory, T violation is equivalent to CP violation.] 



The presence of this unnaturally small number in QCD is known as the strong-CP 
r-| ! problem. 

9 is the sum of two terms, 9qcd and 9qfd- Qqcd is the original value of the 
' angle 9 characterizing the QCD vacuum. 9qfd originates in the chiral rotation 

of the quark fields needed to render the quark masses real and positive. If M p 
and M n are the mass matrices of the up- type (charge 2/3) and down-type (charge 
—1/3) quarks, then 9qfd = arg det(M p M n ). 

There are two general ways of solving the strong-CP problem. In the first 
approach 9 is claimed to have no significance or physical consequences; theories 
with different values of 9 are equivalent, and we may set 9 = without loss of 
generality. This may be so because of the presence in the theory of a Peccei- 
Quinn symmetry [0, but there are also claims that QCD dynamics itself cures 
the strong-CP problem ||. The second path, which I shall follow, tries to find 
some symmetry which naturally leads to the smallness of 9. As 9 is the coefficient 
of a quartic term in the Lagrangian which is both CP- and P-violating, one first 
assumes the Lagrangian — or at least its quartic part — to be either CP- or P- 
symmetric, thereby automatically obtaining 9qcd = 0. CP or P symmetry must 
then be either softly or spontaneously broken. While doing this the problem of 
ensuring the smallness of 9qfd remains. We must obtain quark mass matrices 
which, even though they are complex - - because after all we know that CP 
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is violated in the electroweak interactions — are such that det (M p M n ) is real. 
Afterwards, we must still ensure that the quantum (loop) corrections to the mass 
matrices do not lead to a much too large Qqfd- 

In this letter I put forward a model in which the quartic terms of the La- 
grangian are assumed to enjoy CP symmetry, thereby setting 6qcd = 0. CP 
symmetry is softly broken by the mass terms, which have dimension three, of 
isosinglet vectorlike quarks. The assumption that a symmetry is broken only by 
terms of cubic dimension is technically natural, in the sense that it is not cor- 
rected by infinite renormalizations. Therefore, there will be no non-real divergent 
renormalization of the quartic terms in the Lagrangian. Because of a particular 
discrete symmetry, Oqfd vanishes at tree level. I explicitly show that Qqfd also 
vanishes at one-loop level. My model has some similarity with two previous ones 
|M but is simpler and more economic than any of them. 

The gauge group is SU(2)®U(1) just as in the standard model (SM). The 
scalar sector consists of the usual doublet of the SM (with = ir 2 0*) together 
with a real singlet S. The lepton sector is as in the SM. The quark sector consists 
of the usual three left-handed doublets q^ = (j>l, riL) T , the three right-handed up- 
type quarks Pr, and the three right-handed down-type quarks ur. Besides, I also 
assume the existence of I SU(2)-singlet vectorlike quarks of charge —1/3, with 
left-handed components Nl and right-handed components Nr. The number I is 
arbitrary. 

The model has a discrete symmetry under which both S and the Nr change 
sign while all other fields remain unaffected. As a consequence, the scalar poten- 
tial is 

V = /i!0V + Ax(0V) 2 + ^S 2 + \ 2 S A + A 3 5 2 0V (1) 
and the Yukawa Lagrangian is 

C y = -W C (j>-\n - HZ? ,n - WcStMu +H.\., (2) 

V 

where T and A are 3x3 matrices, T is an I x I matrix. The mass terms —NlMur+ 
H.c, where M is an I x 3 matrix, are relevant too. The vacuum expectation value 
(vev) of 4> is ve ta and the one of S is V. I decompose the scalar multiplets as 

<j) = e ia ( w/5 |, S = V + R. (3) 

v \v+(H + i X )/V2 J 1 ; 

<y9 ± and x are the charged and neutral Goldstone bosons, to be absorbed as the 
longitudinal components of the W ± and Z gauge bosons, respectively. After using 
the vacuum-stability conditions to eliminate fii and \xi in favor of v and V, the 
Higgs potential in Eq. ([I]) yields the mass terms for H and R: 

V = 2X 1 v 2 H 2 + 4X 2 V 2 R 2 + 2V2X 3 vVHR + ... (4) 

It follows that H and R mix to form the two physical scalar particles, which I 
shall denote Si and S%. I denote the mixing angle ip. 
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At tree level 



Q QFD = argdet[{e ia vT){e- ia vA)(VT)} = argdet(rAT) (5) 

does not depend on a. 

Let us denote M u = di&g(m u ,m c ,m t ) the diagonal matrix of the masses of 
the three up-type quarks, and M d = diag(md, m s , m b , ...) the diagonal matrix of 
the masses of the Z + 3 down-type quarks. Without loss of generality, we may 
choose to work in the weak basis in which pi = ul and pn = ur are the physical 
up-type-quark fields, and A = M u /(ve~ ia ). In this basis, we must bi-diagonalize 
the mass matrix of the down-type quarks. This we do by performing the unitary 
transformations (A may be L or R) 

(6) 

in which Xa is a 3 x (Z + 3) and Ya is an Z x (Z + 3) matrix. Indeed, Xl is 
the Cabibbo-Kobayashi-Maskawa (CKM) mixing matrix of the charged gauge 
interactions. I define the (I + 3) x (Z + 3) hermitian matrices Ha = X\Xa- 
As a consequence of the unitarity of the transformation matrices in Eq. @, 
Y\Ya = 1 — Ha, and H\ = Ha- The matrix H L appears in the neutral-current 
interactions of the down- type quarks 0, which are given by 

» . Z ^( s ^pL-l^piH L )d. (7) 





cos 9w \ 3 

The decays of the Z measured at LEP, together with the observed suppression of 
neutral flavor-changing interactions and the near unitarity of the CKM matrix, 
impose strong constraints on the 3x3 part of X^ corresponding to the three 
known charge —1/3 quarks [0]. I shall assume that both M and the vev V are 
very large — in the case of M because it consists of singlet mass terms, in the case 
of V because it is the vev of a singlet field, and it does not break the electroweak 
gauge group. Under that assumption, the 3x3 part of X L is almost unitary, 
as required. Moreover, the mixing angle ip between H and R scales as v/V [see 
Eq. (f|)]. Therefore, in the limit of large V, R is an eigenstate of mass with mass 
~ V. 

The bi-diagonalization condition is 

/ ve ia T 
\ M 

It follows that 

H L M d (l-H R ) = 0, 
H R M d \l-H L ) = 0. (9) 
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The Yukawa interactions of the fields H and R with the down- type quarks are 
given by 

-e ia ^=(nTTn R + H.c.) = -JL(^H L M d d R + H.c), 

—R{NlTN r + H.c.) = -^\d~ L M d {l-H R )d R + H.c], (10) 

where I have used Eqs. (0). 

Up to this point I have not postulated anything about the nature of CP 
violation in this model. I now make the further assumption that the quartic part 
of the Lagrangian conserves CP. Then, Qqcd vanishes, the matrices A, V and T 
are real, and Oqfd = at tree level. CP cannot be spontaneously broken in this 
model, because there is no gauge-invariant phase between the vevs of the scalars 
(the phase a can be gauged away). Therefore, CP must be softly broken by the 
only non-quartic terms in the Lagrangian which may be complex: the matrix M. 
As CP-violating amplitudes are ~ 10 3 times smaller than CP-conserving ones, 
one may naturally assume that the complex phases in the matrix elements of M 
are all ~ 1CT 3 . 

The first model in which soft breaking of CP was used to suppress strong CP 
violation was suggested by Georgi ||. My present model has a clear similarity 
with the one of Bento, Branco, and Parada (BBP) 0]. Both models are of the 
Barr-Nelson || type. In the model of BBP CP breaking is spontaneous and 
there is a complex, instead of real, singlet. There are three neutral scalars, and 
there is CP violation via scalar-pseudoscalar mixing, which leads to strong CP 
violation arising at one-loop level. In contrast, in my model, just as in the one 
of Babu and Mohapatra (BM) M , there are only two neutral scalars, and there 
is no scalar-pseudoscalar mixing, allowing strong CP violation to arise only at 
two-loop level, as I shall now show. The proof is analogous to the one valid in 
the model of BM. 

Strong CP violation might arise at one-loop level via the self-energies of the 
down-type quarks, which give rise to complex mass terms. The relevant diagram 
is the one in Fig. 1. Let us denote the mass term originating from that diagram 
T/f. Ed is an (/ + 3) x (1 + 3) complex matrix. The resulting contribution to 8qfd 
is Im tr(M ( ^" 1 Erf) flQ| . As M d is diagonal, only the diagonal matrix elements of 
Erf are relevant. Now, from Eqs. fllCf ) 

= l ^m k {x{H L M d ) a \H L M d ) kh + y[M d {l-H R )\ ak [M d ^ 

k=l 

+z(H L M d ) ak [M d (l - H R )] kb + z[M d (l - H R )} ak (H L M D ) kb } , (11) 

where the coefficients x, y, and z are real functions, resulting from the loop 
integration, of the squared mass of the k th down-type quark, ml, of the masses 
of Si and S%, and of the mixing angle if). Then, using Eqs. one finds 

1+3 

tr(M^ d ) = J2 m l [(* + *)Hl - (y + z ) r r + vf k ■ (12) 

k=l 
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As Hi and Hp> are hermitian matrices, their diagonal elements are real, and 
Im tr(Mj~ vanishes. 

9qfd therefore arises only at two-loop level. Its magnitude should be sup- 
pressed by at least two factors: a factor (167r 2 ) -2 from the loop integrations, and 
a factor ~ 10 -3 from the small CP-violating phases in M. Other suppression 
factors of the form m q /mw, in which m q is the mass of either the charm or the 



strange quark, may be present too |10| . It should however be noted that in this 
model, such as in the one of BBP, due to the non-diagonal character of Hl, there 
are diagrams contributing to Oqfd at two-loop level with one Z loop and one W 
loop, or with two Z loops. Those diagrams evade the suppressions derived, in 
the context of the SM, in Ref. PH[ . They should however be suppressed by the 
smallness of the non-diagonal matrix elements of Hl. Also, there are diagrams 
with scalar loops contributing to Oqfd at two- loop level, involving in particular 
the mixing matrix Hr [see Eq. fllPD]. Those diagrams will be suppressed if R is 
very heavy (if V is very large), as I advocate. Summing up, it is not clear how 
large may the two- loop contribution to Oqfd be in this model. It is evident any- 
way that, by avoiding scalar-pseudoscalar mixing, my model avoids an important 
source of strong CP violation present in the model of BBP, and in many other 
models. 

I thank Lincoln Wolfenstein for having read a draft of the manuscript and 
commented on it. 
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